TITLE: 
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DESCRIPTION: 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 

AUXILIARY LIBRARY ROUTINE 



M 32 - 29& 



Williams Memory Routine for Linear Programming by 

the Simplex Method -Modified for Variable C. 

Complete Floating Point Program 

This is the same as Library Routine M 15 - 183 "Williams 

Memory Routine for Linear Programming by the Simplex 

Method", except that one of the C.'s may be varied 

o 

by adding an increment AC . . 

All directions are the same as for M 15 - 183 except: 
Prepare the parameter tape as follows: 

005+ 

00F 00 jF 
OOF OOiF 
0019+ 
00F OOtF 
OOkF.OObF 
2^999N 



where 



j = number of activities (Excluding slack and P ) 
i = number of restraints 
t = frequency of intermediate print -out 
k = 

b = for minimization problem 1 for maximization 
problem 



designation of C, to be varied 
k 



In addition, two parameters are needed at the end of the 
data tape. These should be in floating decimal representa- 
tion. The first specifies the increment AG by which C, 
is to be increased for each separate solution. The second 
parameter at the end of the tape specifies the upper limit 

desired for C, . 
k 
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TIMING: 



CAPACITY: 
EXAMPLE: 



The program proceeds by computing an optimal solution 
for the initial value, then it adds £C, to C and 
computes a new optimal solution for C, ' = C v + ^C, • 

K. -K K. 

(The new solution may or may not be different from 
the previous one). It continues in this way until it 
has computed a solution for C, at its upper limit. 

The print-out will be headed by the C, used in comput- 
ing that particular solution. The parameter t controls 
the intermediate print-out, if any, between successive 
final solutions. 

The computational procedure is the same as for M 15 - 183* 
so a comparable amount of time will be required per itera<~ 
tion. However, there will in general be more iterations. 
The time required to change the C. is negligible. 



k(l + j + l) + i j < 682 



Minimize: 



Z = 21^ + 2ky 2 + c^ y^ 

where C^ varies from 2 l/2 to 10 l/2 and Z£L = k 
3 3 

according to the restraints. 

(a) y x -f 2y 2 + y 5 >2 

3y x + 3y 2 + y^ > h 
then (a) can be rewritten: 

(b) yi + 2y 2 + y - y k = 2 



y. >0 





3y 1 + 3y 2 + y^ - 


y 5 


where C, = C =0. 


In vector form this is 


(c) 


12 1-10 

3 3 1 o-l 





y. > 



y 2 




"2 


y 5 




k 


\ 
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but this contains no unit vectors, so (c) is written as 



(d) 



12 1-1 1 o 
3 3 1 0-101 



and 



Z = kfi^ + 24y 2 + c^ + }fy^ + My 

Now if M is set equal to 9 x 10 , for example, then 
it is intuitively obvious that in order for Z to be a 
minimum, y> = y = and thus (c) and (d) would be in 
a sense, "equivalent". 
Then 



P+ 



2 
_k 

-1 




P 2 = P2' 



U 





-1 



and p/r and p_, are the basis vectors 
The coefficient tape is: 



+2+01) 

+3+01J 
+2+0l] 



+3+01 
+1+01 



3 



+1+01 



-i+oD 
+0+00J 



+0+00I 
-1+01/ 



+0+00 ... 


c o 


+21+02") 


C, and C 


+24+02 J 


1 2 


+25+01 


initial C-, 
3 


+0+00") 


c^ = c 5 = 


+0+00J 


+9+O5") 


C 6 = C ? = 9 x 10 


+9+O5J 


+4+01 


AC, = 4 

3 


+105+02 


terminal C^ 



p_ 
9 



4 - 



The parameter tape would be 

005+ 

OOF 005F 
OOF 002F 
0019+ 
OOF OOIF 
003F OOOF 
24999N 

The output has the forms 

+250000000 +01 C^ 

j 

006 +200000000 +01 

006 +400000000 +01 

000 +540000000 +06 

001 +359979000 +06 

002 +449976000 +06 

003 +179997500 +06 

004 -900000000 +05 

005 -900000000 +05 

002 +100000000 +01 

007 +100000000 +01 

000 +900240002 +05 

001 +134991000 +06 

006 -224988000 +06 

003 -449904999 +05 

004 +134988000 +06 

005 -900000000 +05 

002 +666666666 +00 

001 +66666666Q +00 



(5 structural vectors) 
(2 basis vectors) 

(t =D 

to indicate C, varies 
3 

and Z is to be minimized. 
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000 +500000000 +02 
007 -8999^0000 +05 

006 -899970000 +05 

003 +650000000 +01 

004 -299999999 +01 

005 -600000000 +01 

003 +100000000 +01 

001 +100000000 +01 

000 +235000000 +02 

007 -899907500 +05 

006 -900067500 +05 

002 -975000001 +01 
00k +675OOOOOO +01 

005 -925000001 +01 

003 +400000000 +01 

004 +200000000 +01 

000 +100000000 +02 

007 -899975000 +05 

006 -900000000 +05 

002 +65OOOOOOO +02 

001 -135000000 +02 

005 -250000000 +01 

+650000000 +01 c 

003 +400000000 +01 

004 +200000000 +01 

000 +260000000 +02 

007 -899935000 +05 

006 -900000000 +05 

002 -U9999999 +01 

001 -149999999 +01 

005 -650000000 +01 



- 6 - 
+105000000 +02 C 



3 



003 +400000000 +01 

004 +200000000 +01 

000 +420000001 +02 
007 -899895000 +05 

006 -900000000 +05 

002 +750000004 +01 

001 +105000001 +02 

005 -105000000 +02 

003 +100000000 +01 

001 +100000000 +01 

000 +315000000 +02 

007 -899947500 +05 

006 -899947500 +05 

002 +225000002 +01 

004 -525000001 +01 

005 -524999999 +01 

002 +666666667 +00 

001 +666666667 +00 

000 +300000000 +02 

007 -899940000 +05 

006 -899970000 +05 

003 -150000001 +01 

004 -300000000 +01 

005 -600000000 +01 
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00 F 
00 173F 
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00 F 
00 1S5 
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00 3laF 
00 3^1F 
00 S6 
00 S6 



2 


00 F 




00 3^1F 




0021K 





50 12F 




75 6F 


1 


S5 F 




lk 12F 


2 


iA 6f 




^0 7F " 


3 


LO 12F 




LO 6f 


k 


iA 18F 




ho 9F 


5 


iA 12F 




kO 10F 


6 


IA 6f 




ij-0 ]AF 


7 


iA 12F 




kO 15F 


8 


Ik 6f 




^0 11F 



■^LOCATION 


ORDER 


NOTES PAGE 2 M 32 


9 


IX 12F 
1+0 13F 




10 


Lk 12F 
40 8F 




11 


50 12L 
26 999F 




12 


00 F 
00 L 


♦ 




00 F 
26 L 






26lN 


i 


Ik 


LI 6f 
LO 5F 




15 


ko 8l 

kl SF 




16 


F5 SF 
kO ISF 




17 


F5 8L 
kO 8L 




18 


36 9L 
L5 2L 




19 


lA 7L 
40 2L 




20 


26 2L 
00 F 




21 


00 IF 
00 IF 




22 


00 F 
00 F 




23 


50 10L 
26 999F 




2k 


00 F 
00 L 




25 


00 F 
26 L 
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26 
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" 


L5 16F 


27 


42 kh 




50 1L 


28 


26 S4 


; 


OK SN 


: 29 


IK S'6 




8K F 


50 


8S 3F 




15 F 


31 


17 SK 




84 3F 


32 


8S 3F 




12 k-L 


33 


85 3F 




00 S9 


3^ 


OS SS 




8J 69L 


35 


03 3L 




2K S5 


36 


8K f 




8S kF 


37 


21 1SS 




80 kF 


38 


83 l^L 




21 1SS 


39 


8S kF 




8j 71L 


ko 


8J 73L 




23 11L 


kl 


85 kF 




82 l^lL 


k2 


8J 75L 




3K S6 
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8k f 
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44 


83 25L 




31 F 


*5. 


83 25L 




83 25L 


46 


31 34lF 




36 F 


47 


8S 3F 




84 120L 


48 


82 23L 




81 3F 


±9 


83 25L 




81 3F 


50 


8S 120L 




8j 8il 


51 


8 J 73L 




33 17L 


52 


85 120L 




80 121L 


53 


82 65L 




82 65L 


5^ 


4k S6 




45 F 


55 


4s S8 




42 28L 


56 


8j 85L 




5K S6 


57 


8k f 




5S F 


58 


53 31L 




8k IF 


59 


82 33L 




8S F 


60 


6k sn 
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61 


86 20F 
6s SJ 




62 


8J 93L 
62 3^L 




63 


3K SN 
kK S6 




6k 


kl S8 
37 SJ 




65 


kk F 
kS F 




66 


^3 38L 
8J 95L 




67 


32 37L 
5K SN 




68 


55 SJ 
8s F 




69 


8J 97L 
53 ^2L 




70 


8j 99L 
2K S6 




71 


8f IF 
8j 56L 




72 
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89 9F 
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22 k^L 
22 k$L 




7^ 


8F IF 
7K SN 




75 


8j 6ol 
75 SS 




76 


89 9F 
73 ^9L 




77 
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8J 6kL 
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OF F 
50 52L 
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L5 SL 
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82 12F 
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., 



NOTES 



PAGE 8 



M 32 



LOCATION 


ORDER 


NOTES PAGE 9 M 32 


133 


LI 53L 
40 53L 




134 


36 L 
L5 i4f 




135 


L4 13F 
42 103L 


1 
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